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The hyperfine energy levels of alkali metal dimers: 
ground-state polar molecules in electric and magnetic fields 

J. AldegundeU Ben A. Rivington, Piotr S. Zuchowski[!j and Jeremy M. Hutsor|3 
Department of Chemistry, Durham University, South Road, DHl 3LE, United Kingdom 

(Dated: July 11, 2008) 

We investigate the energy levels of heteronuclear alkali metal dimers in levels correlating with 
the lowest rotational level of the ground electronic state, which are important in efforts to produce 
ground-state ultracold molecules. We use density-functional theory to calculate nuclear quadrupole 
and magnetic coupling constants for RbK and RbCs and explore the hyperfine structure in the 
presence of electric and magnetic fields. For nonrotating states, the zero-field splittings are domi- 
nated by the electron-mediated part of the nuclear spin-spin coupling. They are a few kHz for RbK 
isotopologs and a few tens of kHz for RbCs isotopologs. 

PACS numbers: 37.10.Pq, 31.15.aj, 33.15.Pw 
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I. INTRODUCTION 

There is great interest in the formation of ultracold 
molecules and in achieving molecular Bose-Einstein con- 
densation and Fermi degeneracy. Molecules can be 
formed in ultracold atomic gases either by photoassoci- 
ation [1., ^] or by tuning through zero-energy Feshbach 
resonances with magnetic fields [l|, y| . Since alkali metal 
atoms are easier to cool than other species, most work on 
ultracold molecule formation has focussed on alkali metal 
dimers. 

There is particular interest in forming ultracold polar 
molecules. Dipole-dipole interactions are both stronger 
and longer-range than the quadrupole-quadrupole and 
dispersion forces that exist between nonpolar molecules. 
As a result, dipolar quantum gases are predicted to have 
novel properties Q. Ultracold dipolar molecules might 
also be used in quantum information storage and pro- 
cessing [3]. 

Both photoassociation and Feshbach resonance tuning 
form molecules that are initially in highly excited vibra- 
tional states. Quantum gases of such molecules can be 
formed [y, |7J, |8|, but they are long-lived only in very 
specific cases, such as homonuclear fermion dimers in 
the highest vibrational level, tuned to large scattering 
lengths [9]. For other cases the molecules undergo fast 
inelastic collisions that lead to trap loss 10, 11, 12]. Fur- 
thermore, even heteronuclear molecules are essentially 
nonpolar when they are in weakly bound vibrational 
states. Because of this, there is intense current effort 
directed at producing ultracold molecules in their abso- 
lute ground states, for which inelastic losses cannot occur 
and for which heteronuclear molecules have significant 
dipole moments. Very recently, there have been major 
advances in transferring Feshbach molecules to deeply 
bound states by laser-based methods such as stimulated 



Raman adiabatic passage (STIRAP) [13|, [IJ, |15| . Forma- 
tion of quantum gases of ground-state molecules is now 
within reach. 

There has been a considerable amount of work on the 
energy levels of homonuclear alkali metal dimers, espe- 
cially in the near-dissociation states formed by Feshbach 
resonance tuning [1^, Il3, Sa Il3| • However, remarkably 
little is known about the hyperfine structure of the en- 
ergy levels of alkali metal dimers in their lowest rota- 
tional states. The tiny splittings are beyond the resolu- 
tion of most spectroscopic techniques. Nevertheless, an 
understanding of these energy levels is essential in design- 
ing laser-based methods to produce molecules in specific 
states and will be crucial in developing methods to con- 
trol the resulting quantum gases. The purpose of the 
present paper is to investigate the lowest energy levels 
of heteronuclear alkali metal dimers and to explore how 
they behave in electric and magnetic fields. We focus 
here on RbK and RbCs, which are topical for current 
experiments. 



II. THEORY 



Molecular Hamiltonian 
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The Hamiltonian of a diatomic molecule in the pres- 
ence of external magnetic and electric fields can be de- 
composed into six different contributions: the electronic, 
vibrational, rotational, hyperfine. Stark and Zeeman 
terms. By restricting our analysis to ^E molecules in the 
ground electronic state and in a fixed vibrational level, 
the first two terms take a constant value and the rota- 
tional, hyperfine, Stark and Zeeman parts of the Hamil- 
tonian can be written uO, Sl^, S2l 



(1) 



H — Hiot + Hm + Hs + Hz, 



where 



III. EVALUATION OF THE COUPLING 
CONSTANTS 



Hrot 
Hhf 



2 



(2) 



Hz 



i=l 
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Y,C^N■I,+C3Il■T■h+CiIl■h■, (3) 

4 = 1 

-/. • E; (4) 

2 

-g^^jL^N ■ B -'^Qi^T^Ii- B{l-ai). (5) 

i=l 



The three different sources of angular momentum in a ^S 
diatomic molecule are the rotational angular momentum 
N and the spins /i and I2 of nuclei 1 and 2. The rota- 
tional and centrifugal distortion constants of the molecule 
are By and Dy (the centrifugal distortion contribution 
will not be considered in the calculations). The hyper- 
fine Hamiltonian of equation[3]consists of four terms. The 
first is the electric quadrupole interaction with coupling 
constants {eqQ)i and (eqQ)2, where qi is the electric field 
gradient at nucleus i and eQi is its nuclear quadrupole 
moment. The second is the interaction between the nu- 
clear magnetic moments and the magnetic field created 
by the rotation of the molecule, with spin-rotation cou- 
pling constants ci and C2. The two remaining terms rep- 
resent the tensor and scalar interactions between the nu- 
clear dipole moments, with spin-spin coupling constants 
C3 and C4 respectively. The tensor T describes the angle- 
dependence of the direct spin-spin interaction and the 
anisotropic part of the indirect spin-spin interaction [22l | . 

The Stark and Zeeman Hamiltonians, equations [5] and 
[SJ describe the interaction of the molecule with an exter- 
nal electric field E and magnetic field B, where {j, is the 
molecular dipole moment. The Zeeman Hamiltonian con- 
sists of two terms representing the rotational and nuclear 
Zeeman effects. The former arises because the molecular 
rotation produces a magnetic moment grfx^N^ where gr 
is the rotational g-factor of the molecule, which interacts 
with the external magnetic field. The latter arises from 
the interaction of the nuclear magnetic moments gifi^Ii 
with the magnetic field, where gi is the nuclear g-factor 
for nucleus i and li is its nuclear spin. The nuclear shield- 
ing tensor cr^ is approximated here by its isotropic part 
(Ti; terms involving the anisotropy of cr^ are extremely 
small for the states considered here. The diamagnetic 
Zeeman effect is not included in the Hamiltonian as it 
causes level splittings less than 1 Hz for the range of 
magnetic fields considered in this work. 

The nuclear g-factors and quadrupole moments are 
well known [23| . The dipole moments of KRb and RbCs 
have been calculated from relativistic electronic structure 
calculations [H,!!!. 



Nuclear quadrupole coupling constants have been mea- 
sured for several alkali metal dimers as shown in Table 
m However, the only such species for which the mag- 
netic coupling constants have been measured is Na2 j26| | , 
and even there the experiments did not resolve hyper- 
fine splittings for the A^ = state. To the best of our 
knowledge, no experimental data are available for the hy- 
perfine structure of the molecules we consider here, KRb 
and RbCs, in their ground electronic state. We therefore 
carry out electronic structure calculations to estimate 
them. The electric quadrupole coupling constants {eqQ)i 
and {eqQ)2, the nuclear shielding, the spin-rotation con- 
stants Ci and C2 and the spin-spin coupling constants C3 
and C4 are evaluated by density- functional theory (DFT) 
using the ADF package [23, [2^ , which uses Slater func- 
tions and allows the inclusion of relativistic corrections. 
The rotational g-factor (not implemented in the ADF 
code) is evaluated with the DALTON package [28|. 

The objective of the present paper is to explore the 
behaviour of the molecular energy levels in the presence 
of external fields. A detailed discussion of the features 
and effectiveness of the many different methods and basis 
sets available for the calculation of the coupling constants 
is beyond the scope of the work. However, to estimate 
the reliability of the functionals and basis sets employed 
here we compare the coupling constants obtained for a 
group of molecules containing alkali metal atoms with 
experimental results in tables [H |TT1 IIIII and IIVI For sim- 
plicity we have omitted experimental uncertainties and 
vibrational state dependences. It may be seen that the 
calculated coupling constants are generally within 30% of 
the experimental values, except in occasional cases where 
the experimental values are unusually small (such as C4 
for **5Rb35Cl). 

Evaluation of hyperfine coupling constants requires a 
basis set that properly describes the electron density near 
the nuclei. Because of this, we employ all-electron basis 
sets rather than valence basis sets with effective core po- 
tentials. However, for core orbitals of heavy elements 
such as those considered here, relativistic effects can be 
important. In the present work, relativistic corrections 
were included by means of ZORA, the two-component 
zero-order regular approximation |4ll . l42l |43| , including 
spin-orbit coupling as well as scalar effects (which are 
the equivalent of Darwin and mass- velocity terms in the 
Breit-Pauli Hamiltonian) . 

DFT generally performs well for calculations of electric 
quadrupole couplin g co nstants for main-group elements 
Hi m, m, m HM S [U, [H- FoUowingmost of 
these examples, we use the B3LYP functional [5J, [SJ] in 
our calculations with the QZ4P basis set (a quadruple-^ 
all-electron basis set with four polarization functions). 

Shielding tensors were evaluated using the KT2 func- 
tional [55] with the same basis set and relativistic cor- 
rection as for the quadrupole coupling constants. For 



TABLE I: Comparison of electric quadrupole coupling con- 
stants for alkali metals atoms calculated as described in the 
text with experimental values. The units are MHz. 



Molecule 


(egg)C-i- 


{eQqf'^'^ 


Ref. 


23Na2 


-0.456 


-0.459 


[26] 


3^K2 


-0.279 


-0.158 


[30] 


39K19F 


-7.87 


-7.93 


[3^ 


39K^Li 


-0.830 


-1.03 


[32] 


^^K^^Na 


-0.671 


-0.718 


[32] (for K) 


39K23j^a 


-0.216 


0.171" 


[32] (for Na) 


«^Rb2 


-2.283 


-1.1 


[30] 


S5j^jjl9p 


-73.1 


-70.7 


[33] 


ssRb^Sci 


-53.5 


-52.8 


m 


85Rb™Br 


-46.8 


-47.2 


[35] 


85j^jjl27j 


-39.6 


-58.9 


[36] 


«5Rb^Li 


-8.04 


-9.12 


[32] 


133csl9p 


1.30 


1.25 


[37] 


133cs35q 


1.05 


< 1.1" 


[38] 



"Only the absolute value was determined experimentally. 



TABLE II: Comparison between spin-rotation coupling con- 
stants calculated as described in the text and experimentally 
measured. The label 1 refers to the less electronegative atom 
(K, Rb or Cs) and the label 2 to the more electronegative one. 
The units are kHz. 



Molecule c\ 



Calc „Exp Calc „Exp 



-1 



C2 



Ref. 



^^Na2 0.299 0.243 0.299 0.243 [26] 

^^K"F 0.235 0.270 17.5 10.7 [31] 

*^Rb"F 0.598 0.498 16.1 10.6 [33] 

^^Rb^'^Q Q 457 Q 395 q^qq 0394 ^ 

i33Qgi9p ^ Q5 Q gg2 21.9 15.1 [37] 



calculation of shielding tensors of main-group atoms (H, 
C, N, O and F), the performance of this fmictional is 
excellent, and is better f56| than that of more popular 
functional such as BLYPfsl, [13 and B3LYP. 

Two nuclear magnetic moments can interact both di- 
rectly (through space) and indirectly (via the electron 



TABLE HI: Comparison between spin-spin coupling constants 
calculated as described in the text and experimentally mea- 
sured. The units are kHz. 



Molecule C3 



Calc 



Exp Cal 

C3 C4 



^Exp 



Ref. 



^3Na2 0.298 0.303 1.358 1.067 [26] 

^^K^^F 0.470 0.540 0.032 0.030 [31] 

«5Rb"F 0.751 0.797 0.151 0.237 [33] 

^'^Rb^'^Q Q 032 0.033 0.010 0.026 [34] 

i33Qgi9p Q gy5 Q 927 0.471 0.627 [37] 



TABLE IV: Comparison between rotational g-factors calcu- 
lated as described in the text and experimentally measured. 



Molecule g^ 



9: 



Exp 



Ref. 



^^Na2 0.0324 0.0386 [39] 
^^K2 0.0247 0.0212 [39] 

^^Na^^K 0.0253 0.0253 [40] 



85 



distribution) 
action is i22 



Rb2 0.0082 0.0095 [39| 
^Cs2 0.0051 0.0054 [39] 



The coupling constant for the direct inter- 






'), 



(6) 



where R is the internuclear distance. The indirect inter- 
action is represented by a tensor J [22, [53| with isotropic 
part Jiso and anisotropy A J = Jii — Jj_. The coupling 
constants C3 and C4 are related to the direct and indirect 
components by [23, ISSj 



and 



C3 — -Rdd 



C4 = Jis 



AJ 
~3~' 



(7) 



(8) 



In the present work, C3 and C4 were evaluated from equa- 
tions [6] to [8] with {R~^) — Rc^, where Re is the equi- 
librium distance. The components of J were calculated 
using the same methods as for the quadrupole coupling 
constants, except that the PBE 59] functional was used. 
This functional produced results slightly closer to the ex- 
perimental measurements than KT2 for the molecules 
considered in table IIIII (although the differences were 
small). BLYP performed well for all except Na2, for 
which it gave the wrong sign and order of magnitude; 
it also gave qualitatively different results from PBE and 
KT2 for KRb and RbCs. 

ADF does not calculate spin-rotation constants di- 
rectly. However, the spin-rotation constants are given 
approximately by ^60, 61, 62] 



2mcBygi 



((Till — (7ii_) for i = 1,2, 



(9) 



where nip and m, 



are the proton and electron masses. 
By is the rotational constant, gi is the nuclear g- factor 
and cTj II — ai± is the anisotropy of the nuclear shielding 
tensor <tj. Two approximations underlie this expression. 
First, a quadrupole term has been neglected. Secondly, 
it was obtained in the frame of the non-relativistic theory 
developed by Flygare 60]. However, previous studies [63j 
and our own results (see table HI)) suggest that it can be 
applied reliably in the relativistic case. 



TABLE V: Nuclear properties and coupling constants for the different isotopic species of the KRb molecule. 





3»K85Rb 


«SK«^Rb 


«K»5Rb 


«K«^Rb 


«K»5Rb 


«K»^Rb 


Ik 


3/2 


3/2 


4 


4 


3/2 


3/2 


iB.b 


5/2 


3/2 


5/2 


3/2 


5/2 


3/2 


9K 


0.261 


0.261 


-0.324 


-0.324 


0.143 


0.143 


5Rb 


0.541 


1.834 


0.541 


1.834 


0.541 


1.834 


B„/GHz 


1.142 


1.134 


1.123 


1.114 


1.104 


1.096 


(eQ?)K/MHz 


-0.245 


-0.245 


0.306 


0.306 


-0.298 


-0.298 


(eQ?)Rb/MHz 


-3.142 


-1.520 


-3.142 


-1.520 


-3.142 


-1.520 


crK(ppni) 


1321 


1321 


1321 


1321 


1321 


1321 


crRb(ppm) 


3469 


3469 


3469 


3469 


3469 


3469 


ck/Hz 


19.9 


19.8 


-24.2 


-24.1 


10.5 


10.4 


CRb/Hz 


127.0 


427.5 


124.8 


420.1 


122.8 


413.1 


C3/Hz 


11.5 


38.9 


-14.2 


-48.2 


6.3 


21.3 


C4/Hz 


482.5 


1635.7 


-599.0 


-2030.4 


264.3 


896.2 


5r 


0.0144 


0.0142 


0.0141 


0.0140 


0.0139 


0.0138 


m/d 


0.76 


0.76 


0.76 


0.76 


0.76 


0.76 



Lastly, the rotational g-factors were evaluated with the 
DALTON program using the KT2 functional and the all- 
electron basis sets of Huzinaga and coworkers [GJ, [6a |. 
Again, the choice of the functional is based on its reli- 
ability for this molecular property |66| . No relativistic 
corrections were included in this case. Previous calcu- 
lations [63 for hydrogen halides and noble gas hydride 
cations including atoms as heavy as I and Xe suggest that 
relativistic corrections are relatively small for rotational 
g-factors (less than 5% of the non- relativistic value). 

The coupling constants obtained for KRb and RbCs 
are given in tables |V] and IVII All the calculations were 
carried out at the equilibrium geometries. Re = 4.07 A 
for KRb \6§\ and R^ = 4.37 A for RbCs ^69^. This ne- 
glects small corrections due to vibrational averaging even 
for V — 0, but nevertheless gives results that are quali- 
tatively valid for any low-lying vibrational state. ADF 
generally gives coupling constants for only one isotopic 
species, but the others may be obtained by simple scaling. 
The nuclear quadrupole coupling constants scale with the 
nuclear quadrupoles Qi , the spin-spin coupling constants 
with the product of nuclear g-factors gtgj , and the spin- 
rotation coupling constant with the product of gi and the 
rotational constant By. The rotational g- factor scales in 
a more complicated way that depends on By and the shift 
of the center of mass [TOl • 



TABLE VL Nuclear properties and coupling constants for the 
different isotopic species of the RbCs molecule. 





85Rb"3Cs 


^^Rb^^^Cs 


^Rb 


5/2 


3/2 


/Cs 


7/2 


7/2 


ffRb 


0.541 


1.834 


gcs 


0.738 


0.738 


B„/GHz 


0.511 


0.504 


(eQg)Rb/MHz 


-1.803 


-0.872 


(eQq)cs/MHz 


0.051 


0.051 


crRb(ppm) 


3531 


3531 


o-cs(ppm) 


6367 


6367 


CRb/Hz 


29.4 


98.4 


ccs/Hz 


196.8 


194.1 


cs/Hz 


56.8 


192.4 


C4/HZ 


5116.6 


17345.4 


ffr 


0.0063 


0.0062 


M/D 


1.25 


1.25 



basis sets. 



\I1M1I2M2NMN) 

\{hh)IMiNMN) 

\{IiI2)INFMf) 



(uncoupled basis); (10) 

(spin-coupled basis); (11) 
(fully coupled basis). (12) 



IV. HYPERFINE ENERGY LEVELS 



We calculate the hyperfine levels by diagonalizing the 
complete Hamiltonian of equations [D to [5] in a basis set of 
angular momentum functions. We employ three different 



Here / and F are quantum numbers for the total nuclear 
spin and total angular momentum and Mj and Mp rep- 
resent their projections onto the Z axis defined by the 
external field. We consider here only cases in which only 
one field, electric or magnetic, in present. The matrix ele- 
ments corresponding to the different terms of the Hamil- 



tonian in each of the basis sets are calculated through 
standard angular momentum techniques [7l| . 

The use of three basis sets rather than one helps in as- 
signing quantum numbers to the energy levels. Although 
the Hamiltonian matrix is not diagonal in any of the ba- 
sis sets employed, it is usually closer to diagonal for one 
basis than for the others. When one coefficient of an 
eigenvector is much larger than the others, it is possi- 
ble to assign approximate quantum numbers to the state 
concerned. However, different basis sets achieve this in 
different field regimes. 



Zeeman splitting for rotational ground-state 
molecules (TV = 0) 



Figure [T] shows the Zeeman splittings for energy levels 
of ^^K^^Rb with iV = 0. The splittings are dominated by 
the scalar nuclear spin-spin interaction and the nuclear 
Zeeman effect, which are the only terms in the Hamilto- 
nian with matrix elements diagonal in N for A^ = 0. It 
should be noted that the scalar spin-spin coupling is en- 
tirely mediated by the electron distribution, and has no 
contribution from the direct dipolar interaction. In the 
absence of external fields, the energy levels are split into 
groups labeled by the total nuclear spin /. For small mag- 
netic fields B, I remains a nearly good quantum number 
and the levels split according to the value of its projec- 
tion Mj (which in this case coincides with the projection 
of the total angular momentum, which is always a good 
quantum number). Energy levels corresponding to the 
same value of Mj display avoided crossings as a function 
of the field as shown in figure [H For fields well above 
the crossings (which are at 2 to 10 G in this case), / 
is destroyed and the good quantum numbers are M^b 
and Mk. Since both nuclear g- factors are positive for 
^^K^^Rb, states where both projections are positive are 
high-field-seeking and those where both are negative are 
low-field-seeking. 

Although the splittings at low fields are dominated by 
the scalar spin-spin coupling, there are several terms in 
the Hamiltonian that are off-diagonal in A^ . The energies 
are therefore obtained by diagonalizing a full matrix that 
includes enough rotational levels for convergence. For 
the Zeeman effect, the only off-diagonal terms involving 
A^ = are the electric quadrupole coupling and the tensor 
spin-spin coupling, both of which are small. Convergence 
for A^ = is achieved with A^max = 2 and the splittings 
obtained differ from those calculated with only A^ = 
by less than 1%. For the Stark effect, however, the Stark 
term itself mixes A^ = states with A^ > 0. Terms off- 
diagonal in A^ are then very important and much larger 
basis sets are needed. 

The scalar spin-spin interaction for A^ = is diagonal 
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FIG. 1: Zeeman levels for 39K*5Rb(« = 0, A = 0). 
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FIG. 2: Zeeman splitting and avoided crossings (indicated 
with arrows) for the Mi = -1 levels of ^^K**^Rb(t; = 0, A = 
0). 



in the spin-coupled and fully coupled basis sets, 

(AT = 0(/i/2)/Af/|c4/i -HN ^ 0{hl2)IMi) = 
{N = 0{hl2)IFMF\ciIi ■ hlN = Q{Iih)IFMF) 

ic4[/(/ + 1) - h{h + 1) - h{h + 1)]. (13) 

The nuclear Zeeman Hamiltonian is diagonal in the un- 
coupled basis set, with nonzero elements given by 

- [5iA/i(l - (Ti) + .92Af2(l - <J2)]mB. (14) 

The splitting pattern is therefore determined by the al- 
lowed values of the total nuclear spin quantum number / 
and by the magnitudes and signs of the scalar spin-spin 
coupling constant C4 and the rotational g-factors. The 
nuclear shielding constants ai are only a few parts per 
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Magnetic Field (G) 

FIG. 3: Zeeman levels for **5Rb"3Cs(ii = 0, iV = 0). 



thousand. For large values of the magnetic field, where 
the nuclear Zeeman effect is the dominant term in the 
Hamiltonian, the magnetic moment (gradient of the en- 
ergy with respect to B) is close to —{giMi + g2M2)fJ.N- 

The Zeeman splittings for ^'^Rb^'^^ Cs are shown in 
figure [3l They are qualitatively similar to those for 
•^^K^^Rb, except that the range of / is different and the 
spin-spin coupling constant C4 is significantly larger. Be- 
cause of this, I remains a good quantum number up to 
significantly higher magnetic fields. At high fields, once 
the magnitude of the scalar spin-spin interaction can be 
neglected compared to the Zeeman effect, MRb and Mcs 
become good quantum numbers. 

The splitting patterns for other KRb and RbCs iso- 
topologs are qualitatively similar to those discussed 
above and the corresponding figures are available as 
supplementary online material. The spin-spin coupling 
constant and the potassium g-factor are negative for 
*°K*^Rb and ^OxS^Rb. The sign of a determines 
whether the lowest zero-field energy corresponds to the 
highest or lowest value of /. In general the fields where 
the avoided crossings occur and above which Mi and M2 
become good quantum numbers scale with {04/ {gi — g2)\- 
When gi and 52 are equal, as in homonuclear dimers, 
there are no avoided crossings for N = and the / quan- 
tum number is conserved even at high fields. 



B. Stark splitting for rotational ground-state 
molecules (N = 0) 

The Stark effect for levels of ^^K^^Rb correlating with 
A^ = is shown in figure [J] to [51 Corresponding fig- 
ures for the remaining isotopologs of KRb and RbCs are 
available as additional online material. The Stark ef- 
fect is quadratic at low fields but becomes linear at high 
fields, as is usual for diatomic molecules in E states [72l |. 
This arises from mixing between different rotational lev- 
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FIG. 4: Stark effect on energy levels of K Rb correlating 
with {v — 0,N = 0) for electric fields up to 50 kV/cm. 
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FIG. 5: Stark splitting for energy levels of '^^K^^Rb correlat- 
ing with {v = 0,N = 0) for electric fields up to 1 kV/cm. 
The levels are shown relative to their field- dependent average 
energy. 



els: while in the Zeeman case this mixing is very weak 
and is exclusively due to hyperfine terms, in the Stark 
case it is strong and is caused directly by the electric 
field. At low fields the mixing is weak and can be treated 
by second-order perturbation theory, giving rise to a 
quadratic Stark effect. However, as the field increases 
the mixing becomes increasingly important: the A^ = 1 
basis functions contribute around 25% at 10 kV/cm and 
40% at 20 kV/cm. Eventually the molecule becomes fully 
oriented by the field and the linear Stark effect overcomes 
the quadratic effect. The mixing also has numerical con- 
sequences as the number of rotational levels required for 
convergence increases with field: for example, calcula- 
tions at 50 kV/cm require A^max = 6. 

The magnitude of the Stark shift in figure 2] obscures 
the splittings between hyperfine levels. Figure[5]therefore 
shows the levels correlating with A^ = relative to their 
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average energy, for fields up to 1 kV/cm. As expected, 
each zero-field level splits into I + 1 components labeled 
by the different possible values of |M/|. For |M/| > the 
levels exist in degenerate pairs corresponding to chang- 
ing the sign of Mi and M2. However, changing the sign 
of one of Ml and M2 produces a different state with a 
different value of |A//|. For AIj — there is an extra sym- 
metry corresponding to reflection in a plane containing 
the electric field vector. 

At higher field, as shown in figure [SI the projections of 
the individual nuclear spins become well-defined as well 
as their sum. At sufficiently large fields the splittings ap- 
proach a limiting value as the molecules become strongly 



oriented along the field direction. In this limit the split- 
tings are mostly determined by the nuclear quadrupole 
coupling constants, with relatively small contributions 
from the magnetic hypcrfinc terms. 



V. CONCLUSION 

We have investigated the hyperfine level splittings ex- 
pected for alkali metal dimers in their rotational ground 
state in the presence of electric and magnetic fields. We 
have carried out density-functional calculations of the 
electronic structure of RbK and RbCs at the equilibrium 
geometry of the ground ^S state and evaluated all the 
hyperfine coupling constants necessary to calculate en- 
ergy level patterns. For nonrotating states, the zero-field 
splittings between hyperfine states range from a few kHz 
for isotopologs of KRb to a few tens of kHz for isotopologs 
of RbCs. They are dominated by the electron-mediated 
contribution to the nuclear spin-spin coupling. The re- 
sults will be valuable in designing laser-based schemes 
to produce ultracold molecules in their absolute ground 
states in applied fields. 
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